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NYO-7697 

BASIC  NOTIONS  OF  RELATIVISTIC  HYDROMAGNETICS 

Paul  Relchel 

Part  I  -  Formulation 

1.   Introduction 

The  piorpose  of  Part  I  is  to  formulate  the  laws  of  hydro- 
magnetics  in  the  covariant  four  dimensional  form  of  special 

[2] 
relativity,   Taub    has  formulated  the  corresponding  relativis- 

tic  hydrodynamic  equations  for  the  non-electromagnetic  case,,   On 
the  other  hand  de  Hoffman  and  Teller    have  set  up  shock  con- 
servation conditions  (jump  conditions)  for  the  relativistic  treat- 
ment of  hydromagnetic  shocks.   These  conservation  conditions  are 
not  adequate  to  permit  of  actual  solutions  for  the  hydromagnetic 
field  structure,  but  determine  only  shock  transition  values. 
The  present  formulation  extends  the  work  of  de  Hoffman  and  Teller 
in  that  it  consists  of  differential  equations  determining  the 
hydromagnetic  field  structure  (when  boundary  and  initial  condi- 
tions are  given).   Thus  the  conditions  of  de  Hoffman  and  Teller 
can  be  derived  from  the  present  equations. 
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2,   Method  of  Obtaining:  Formulation  -  A  Search  for  the  Simplest 
Covarlant  Expressions 

The  derivation  procediire  of  a  tensor  equation  (I.e.  Lorentz 
invariant  in  a  (+++-)  space)  is  tantamount  to  first  selecting  or 
constructing  "promising"  four -dimensional  tensors  and  then  examin- 
ing these  to  see  If  the  laws  of  physics  may  be  expressed  in  terms 
of  equations  involving  them.   An  outstanding  guiding  rule  in  this 
procediire  (for  the  case  of  special  relativity  only)  is:   if  we 
can  express  merely  the  non-relativistic  laws  of  physics  (i.e. 
equations  describing  low  velocity  and  low  energy  effects)  in 
terms  of  tensors,  these  equations  often  constitute  automatically 
the  desired  invariant  forms. 

More  specifically,  tensor  equations  may  be  found  largely  by 
two  techniques: 

(a)  A  classical  equation  may  frequently  be  enhanced  into  a 
tensor  equation  essentially  by  "adding  an  appropriate  extra 
dimension"  to  the  three  dimensional  "vectors"  or  matrices  appear- 
ing in  the  original  equation;  and  then  checking  via  a  Lorentz 
transformation  that  such  a  generalization  actually  constitutes  a 
tensor  equation.   This  usually  means  the  four  dimensional  vectors, 
matrices,  etc.,  entering  the  equation  are  Individually  tensors. 

(b)  In  the  procedure  of  (a),  roughly  speaking,  we  selected 
tensors  to  fit  an  equation.   Frequently  the  reverse  is  done: 
starting  with  some  tensor  an  equation  is  "built  around  it."   It 
is  worth  noticing  that  unlike  in  the  former  case,  the  present 
tensor  need  not  necessarily  have  been  foTxnd  by  conjecture,  but 
may  have  been  directly  constructed  via  operations  on  previously 
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known  tensors. 

While  physical  and  mathematical  considerations  are  of  great 
value  in  suggesting  possible  covariant  expressions,  it  is  clear 
from  the  foregoing  that  the  actual  derivation  consists  of  little 
more  than  finally  recognizing  that  a  trial  expression  does  in 
fact  transform  covariantly,  and  does  have  significant  physical 
content.   Thus  the  formulation  will  merely  be  stated  and  examined 
with  no  attempt  at  a  systematic  derivation. 
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3,   Physical  Assumptions  (A  Relatlvistic  Extension  of  Lundquist's 
Equations) 
In  order  to  permit  reasonable  mathematical  handling  the 
equations  to  be  formulated  v-rill  describe  an  idealization  of  a 
somewhat  specialized  system  frequently  occurring  in  nature. 
The  following  specific  limitations  will  apply: 

(a)  The-  medixAm  in  question  will  consist  of  several  ionic 
fluid  components.   We  will  consider  only  the  case  where  the 
relative  velocities  (diffusion  velocities)  at  any  given  point 
between  these  fluid  components  is  small  compared  to  the  speed  of 
light.'"'   (However  the  relative  velocities  between  different  points 
in  the  fluid  need  not  be  small  compared  to  the  speed  of  light.) 
This  assumption  is  necessary  in  order  that  the  system  may  be 
described  and  solved  in  terms  of  the  total  fluid  properties  with- 
out regard  to  the  individual  component  fluids.   Roughly  speaking, 
we  may  say  the  system  has  been  reduced  to  a  single  effective  fluid 
with  a  "massless  electric  current  flow  tacked  on." 

(b)  We  assume  the  flow  to  be  adiabatic  (e.g.  zero  heat  con- 
ductivity) . 

(c)  Zero  viscosity  is  assumed. 

(d)  We  ass-ume  what  is  frequently  termed  "Ohm's  law  for  a 

perfect  conductor":   E  +  u  x  B  =  0,   This  relation  may  fail  to 

describe  actual  dissipationless  media  in  particular  when  space  or 

time  derivatives  of  the  hydromagnetlc  field  variables  are  too  large. 

"For  example,  if  in  a  two-component  fluid  the  components  are  of 
say  equal  mass  densities  and  travel  at  large  eaual  and  opposite 
velocities  through  the  same  region  simultaneously,  then  the 
velocity  of  the  total  fluid  is  zero  and   though  each  fluid  compo- 
nent is  undergoing  a  Lorentz  contraction,  the  factor 
/f      A) — ^  ' 

i-u  /c   will  be  iinity.   This  can  be  seen  to  invalidate 
Equation  (6). 
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Nevertheless  it  has  the  advantages  of  applying  In  many  cases, 
of  simplicity,  and  of  Lorentz  invariance  of  form^ 

(e)  The  pressure  (relative  stress)  is  assumed  to  be  isotropic, 

(f)  We  assume  the  dielectric  constant  and  permeability  of 
the  medium  to  approximate  that  of  empty  space,   (Hence  |j,  X  =  c°^o) 
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The  3-space  components  of  fluid  velocity:  u  ,  u  ,  u  , 

X   y   z 


I4.,   Symbols 

All  definitions  refer  to  the  flat  space  time  of  special 

relativity. 

11    22    33 
The  metric  tensor:     g   =  g   =  g   =  1 

811  =  822  "  833  =  1 

All  off -diagonal  elements  are  zero. 

The  3-space  components  of  fluic 

Conserved  rest  mass  per  xxnit  3-space  proper  volume  =  Pq» 
The  conserved  rest  mass  is  defined  as  the  total  rest  mass  less 
the  part  we  regard  as  "internal  energy."   (The  meaning  of  a  "true 
zero  level"  of  internal  energy  is  not  important  in  this  connection. 
It  is  only  necessary  that  we  designate  some  standard  thermodynamic 
state  as  corresponding  to  zero  internal  energy.   The  "conserved 
rest  mass"  is  in  fact  conserved  only  if  we  do  not  consider 
creation  or  annihilation  of  fundamental  particles  -  e.g.  pair 
production. ) 

The  "four  space  velocity"  vector  (as  it  is  often  called), 
K  or  K^:   for  a  =  1,2,3:   K°^  =  K  =    ^ 


/^ 


K^^ 


1-^ 


2 
1 

-; 

c' 
2 


-  9  - 


Proper  internal  energy  per  unit  conserved  rest  mass  =  e 
Scalar  invariant  pressiare  =  p. 


o 


Proper  enthalpy  per  unit  conserved  rest  mass  =  i  =  e  +  ^, 
■^  *  -   o  -   o   p^ 

o 


Electromagnetic  tensor:   B  ,,  =  -Bj,  . 
as  the  following  array: 

r 


This  may  be  written 


0 


-B 


B 


B 


-B. 


X 


-B 
3 

0 


E., 


E_ 


E 


-E    -E     -E     0 

If  this  is  changed  to  contravar iant  form  the  result  is  as  the 
above  array  except  that  the  elements  forming  the  fourth  rovj 
and  column  are  each  multiplied  by  -c~     /  -  U-  -    -   . 


B 


h^ 


(b'^^  =  -c""B  ,  and 


=  -c"^B,  .  ) 

Charge  density  (seen  in  any  given  Lorentz  frame)  =  q, 
Cxorrent  fo\ar  vector:   I  .   Fourth  contravar  iant  component  = 


q/cS 


Permeability  =  jx  and  dielectric  constant  =  X.  • 

9 


Coordinate  derivative  operator  S  V  - 


9x, 


a 


10 


5o   Covarlant  Formulation 


(1)       K^K^  =  -c^  , 

This  takes  the  form  of  an  Identity  If  we  regard  the  K' s  as 
defined  in  terras  of  u  ,  u  ,  u  : 

2    .      v2     ,        2     ,       ,2 


^Fij    f^l    \f^, 


Thus,  in  the  three  dimensional  formulation  where  the  u  ,  u  ,  u 

X    y    z 

variables  are  included,  (1)  will  be  an  identity;  while  in  the 
equivalent  four  dimensional  formulation  where  these  are  not  in- 
cluded, (1)  will  be  a  constraint, 

(2)  (REFERS  TO  A  GIVEN  ELEMENT  OF  FLUID):     p  =  p(Pq) 

This  general  expression  for  the  equation  of  state  indicates 
that  only  one  thermodynamic  variable  is  needed  to  determine  the 
state,  in  accordance  with  our  assumption  that  the  flow  is  adiabatlc, 
While  this  form  of  equation  applies  to  each  fluid  element,  we 
have  not  assumed  that  all  elements  obey  the  same  particular 
equation  of  state,   (i.e.   The  formulation  is  not  restricted  to 
a  system  in  which  the  entropy  per  unit  conserved  rest  mass  is 
space  independent.) 

(3)  B^pRp   =   0  . 

This  relation  is  the  covarlant  form  of  what  is  known  as 
"Ohm's  law  for  a  perfect  conductor",  and  states  that  an  observer 
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stationary  with  respect  to  the  fluid  sees  zero  electric  field. 
Equation  (3)  is  completely  equivalent  to  the  following  three 
dimensional  expression: 

(3A)       E+uxB=0 
(3B)       E  •  u  =  0  . 

Thus  (3B)  merely  states  that  E  is  perpendicular  to  u,  a 
fact  already  obvious  from  ( 3A ) . 

If  (3)  be  premultiplied  by  K  ,  then  because  of  the  skew 
symmetry  of  B  „,  the  result  is  identically  zero: 

(3c)     k^'b^-kP  =  0 

or   K^(B^pKP)  +  k2(B2pKP)  +  K^(B3pKP)  +  K^Ib^^rP)  =  0. 
Since  the  four  brackets  are  respectively  the  four  expressions 
represented  by  (3)  we  see  that  any  one  of  these  foiir  expressions 
is  linearly  dependent  on  the  other  three.   This  then  constitutes 
an  alternate  proof  that  {3B)  is  merely  a  consequence  of  ( 3A ) . 

(I;)       V  B„   +  V7„B    +  y  B  „  =  0. 
^        Va  Py    v(3  ya         Vy  ap 

This  is  a  tensor  eouation,  for  cyclic  differentiation  of  a  skew 

symmetric  tensor  yields  a  skew  symmetric  tensor  one  rank  higher 

(in  the  flat  space  time  of  special  relativity).    Equation  (Ij.)  is 

equivalent  to  two  of  the  Maxwell  equations: 

(Ua)       curl  E  +  II  =   0 

ot 

(1|B)       div  B  =  0, 

Equation  (I4E)  is  linearly  independent  of  the  three  equations 
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comprising  (Ua),  but  In  a  rather  trivial  way.   Thus  ( I4A )  Implies 

■ps 

not  (UB),  but  the  weaker  statement  that  -5-7  dlv  B  =  0.   However, 

0 1 

if  we  assume  as  a  supplementary  condition  that  at  any  one  time 

dlv  B  =  0  at  a  point  then  it  will  always  be  so  at  that  point. 

that  is  (I4B)  follows  from  ( I4A )  combined  with  the  physically 

reasonable  assumption  that  div  B  vanished  at  some  time  or  other 

in  each  part  of  the  region  of  interest. 

The  dependence  between  the  "weakened"  ( UB ) ,  namely  that 

—  div  B  =  0,  and  (Ua)   may  also  be  seen  by  cyclically  differen- 
dt 

tiating  {l\.) : 

ikC)  V,B^^)  -  V^(  V,Bp,-.  VpB^^.  VgB^p)  - 

y^  (  V  B„  +  v„B   +  V  B  „ )  =  0  . 

This  is  seen  to  be  an  Identity  because  of  the  skew  symmetry  of 
B„Q«   Taking  5  =  li,  the  fourth  term  becomes  ^r—  div  B,  which  is 
seen  to  be  linearly  dependent  on  space  derivatives  of  the  three 
equations  comprising  (IiA), 

(5)      v^-^    VpB^P  =  i^  . 

This  is  equivalent  to  two  of  the  Maxwell  equations: 

(5a)      curl  H  =  1^  +  I 

dt 

(5B)      div  D  =  q   {=  c^I^)  . 

Equation  (5B)  and  the  three  equations  comprising  ( ^A )  are 
linearly  independent.   Eouatlon  (5)  will  not  enter  explicitly 
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into  the  formulation,  but  will  be  eliminated  by  replacing  I 
wherever  it  appears  by  ^l"  VrB°'^.   It  is  however  interesting  to 
examine  (5)  briefly: 

For  the  I*^  to  form  a  non-zero  four  vector  it  is  readily 
seen  that  we  cannot  have  q  =  0  in  all  reference  frames.   In  other 
words,  except  for  the  case  of  I  =  0,  it  is  impossible  for  a  fluid 
to  be  neutral  in  all  reference  frames.   To  understand  this 
physically,  we  may  consider  a  two  component  fluid  in  which  the 
positive  part  Is  moving  at  some  given  velocity  with  respect  to 
the  negative  part.   Then  an  observer  stationary  with  respect  to 
the  negative  part  sees  the  positive  part  moving  past  him  and 
therefore  of  increased  charge  density  due  to  a  Lorentz  contrac- 
tion.  On  the  other  hand  an  observer  stationary  v/ith  respect  to 
the  positive  charge  does  not  see  it  londer  a  Lorentz  contraction, 
but  instead  sees  the  negative  charge  "compressed".   The  fact  that 
different  observers  may  disagree  on  the  relative  number  of  posi- 
tive and  negative  particles  in  a  given  volTime  element  of  the 
medium  does  not  constitute  an  inconsistency,  but  is  in  fact  a 
consequence  of  the  different  meanings  of  simultaneity  for-  the 
different  observers.   Clearly  even  if  the  diffusion  velocity  is 
small  compared  to  that  of  light,  a  medi'urn  may  display  opposite 
net  charge  signs  in  two  slightly  different  Lorentz  frames. 

Equation  (5)  (or  what  is  the  same  thing,  ( 5A )  and  ( 5B ) ) 
implies  the  conservation  of  charge;  for  by  taking  the  divergence 
of  (5a)  and  the  time  derivative  of  (5B),  the  E  variable  may  be 
eliminated,  yielding  the  relation: 
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61    dl    dl^ 
dt       dx    oy    dz 

(6)       ^a^PoK""^  =  °  • 

Eauation  (6)  states  the  conservation  of  the  Quantity  we  have 
accordingly  called  "the  conserved  rest  massj'   This  relation  has 
been  used  by  Taub  in  the  corresponding  non-electromagnetic 
formulation  (reference  [2]),   In  terras  of  the  fluid  velocity  u 
it  becomes: 

Vt  +  V  .    ^    =  0. 

Remembering  that  p   is  not  the  conserved  rest  mass  per  unit 

volume  seen  by  an  observer  whom  the  fluid  is  passing  at  speed  u, 

but  is  defined  as  the  conserved  rest  mass  density  in  the  proper 

frame  (I.e.  stationary  with  respect  to  the  fluid),  we  recognize 

2  -1/2 
the  necessity  of  the  factors  (1-  — ^)      in  { 6a ) .   This  conserva- 

c 
tion  law  must  be  Qualified  as  follows: 

(a)  There  are  notable  exceptions  to  this  conservation  law. 
For  example  if  an  electron  and  positron  annihilate  each  other  by 
forming  electromagnetic  radiation,  the  original  rest  mass  simply 
disappears,  since  photons  have  zero  rest  mass.   Similar  exceptions 
are  to  be  found  in  the  various  nuclear  reactions,  etc.   There 

are  actually  two  conservation  laws  instead  of  one  (to  the  same 
approximation):   energy  and  the  number  of  particles, 

(b)  In  order  that  a  conservation  law  apply,  the  quantity 
to  be  designated  as  "conserved  rest  mass"  must  by  definition  be 

set  aside  from  the  internal  energy.   Thus  electromagnetic  radiation 
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which  has  zero  rest  mass,  can  supply  internal  energy  to  a 
material  (by  absorption  or  adiabatically  by  say  a  "radiation 
pressiire  compression");  so  that  including  the  internal  energy 
in  the  definition  of  conserved  rest  mass  would  in  fact  violate 
the  desired  conservation  property.   It  need  not  be  supposed  how- 
ever that  there  is  some  inherent  "true  zero  level"  of  internal 
energy,  so  that  the  total  proper  mass  may  be  uniquely  subdivided 
into  a  conserved  rest  mass  part  and  an  internal  energy  part.   It 
is  sufficient  merely  to  regard  the  conserved  rest  mass  as  the 
total  proper  mass  of  some  arbitrarily  chosen  "standard  state", 
and  then  to  interpret  the  internal  energy  of  any  other  state  as 
either  positive  or  negative  depending  upon  whether  the  total 
proper  mass  Is  respectively  more  or  less  than  in  the  standard 
state.   (Negative  internal  energy  is  a  reasonable  concept:   For 
example  at  0°  Kelvin  the  negative  lattice  binding  energy  of  a 
material  far  outweights  the  slight  positive  zero-point  vibrational 
energy. ) 

Since  (6)  permits  of  the  above  very  arbitrary  definition  of 
conserved  rest  mass  -  as  some  constant  part  of  the  variable  total 
proper  mass  -  it  would  appear  that  the  equation  is  too  "artificial" 
to  be  of  much  virtue.   In  other  words  it  may  be  asked  why  do  we 
not  simply  reduce  the  number  of  variables  and  eouations  by  one 
by  removing  the  conserved  rest  mass  and  enthalpy  variables  and 
attempting  to  employ  instead  a  single  total  proper  mass  variable, 
and  dropping  (6)?   Actually  the  importance  of  the  concept  of 
conserved  rest  mass  arises  because  In  relativistic  hydrodynamics 
it  appears  to  be  perhaps  the  unique  quantity  through  which  we  may 
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quantitatively  designate  a  discrete  element  of  fluid  (unless 
we  borrow  microscopic  concepts  to  designate  it  -  e.g.  ions  enter- 
ing and  leaving  the  element).   Thus  a  discrete  element  of  fluid 
is  defined  as  the  fluid  within  a  closed  siirface  which  moves  so 
that  there  is  no  net  flow  of  conserved  rest  mass"  across  it  at 
any  point.   From  this  definition  (and  (6))  we  see  that  the  con- 
served rest  mass  of  an  element  of  fluid  is  a  constant  of  the  motion. 
Hence  the  mathematical  usefulness  of  the  concept  of  conserved  rest 
mass  arises  because  it  enables  us  to  use  a  convenient  type  of 
Lagrangian  derivative  formulation;  for  without  the  conserved  rest 

mass  quantity  a  Lagrangian  formulation  would  require  "semi- 

ENTHALPY 


quantitative"  variables  like  (■ 


) ,  but  with  the 


■GIVEN  FLUID  ELEMENT 
help  of  the  conserved  rest  mass  quantity  this  variable  assumes  the 

ENTHALPY 


quantitatively  meaningful  form  ( 


UNIT  CONSERVED  REST  MASS 


). 


(7) 


p^kP  y  (i+c-2i^)K^.  V-^^^"  \v/^  =  ° 


This  is  equivalent  to  the  following  two  equations  (using 
the  fact  that  we  have  assiomed  [i  ?C  =  c   ): 


(7A) 


lo_  A. 

— o  dt 


fl- 


u 
c 


(1+  -^)—^ 
c"^  /    2 


1- 


u 


Vp+bx  Iv  xh-V^d| 


-E(V-D)  =  0 


(7B) 


'^o    d 

TT  dt 


/1-H, 


+  1^  +  E  o  (VxH  -  V^d|  =  0, 


"Microscopically  this   means   no   net   flow  of  particles, 
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Eauatlon  (7A)  is  the  equation  of  motion  corresponding  to 
Newton's  second  law.   Equation  (7B)  describes  (but  not  very 
obviously)  the  rate  of  change  of  the  total  (non-electromagnetic) 
mass  energy  of  an  element  of  fluid  due  to  a  combination  of  mechani- 
cal and  electromagnetic  effects.   In  the  proper  frame  we  have 
from  (3)  that  E  =  0,  so  that  (7B)  reduces  In  that  frame  to 


o  dt   o    dt 


or 


(7C)     dl^  =  |E  . 

^O 

This  relation,  which  Is  a  consequence  of  (3)  and  (7B), 

states  that  the  entropy  of  a  given  element  of  fluid  (appropriately 

defined)  remains  unchanged.   (Using  the  definition  1   =  e  +  pp~  . 

o  -   o   ^^o  * 

(7C)  can  readily  be  shown  to  describe  an  element  of  fluid  when  the 

only  Internal  energy  changes  are  due  to  compression  -  not  heat 

conduction  or  diffusion.)   This  will  render  the  equation  of  state 

an  "adiabatic"  one,  ^ 

Equation  (73)  is  not  a  consequence  of  ( 7A )  alone,  but  will 

now  be  shovjn  to  be  a  conseauence  of  (1),  (3),  ( 7A )  and  ( 7C ) .   If 

(7)  is  premultiplied  by  K  we  get  a  linear  combination  of  the 

foup  equations  represented  by  (7).   Performing  the  indicated 

differentiations  and  putting  di^  =  p"  dp  ,  this  linear  combination 

is 

p^(l^c-2i^)KPK-  VpVK-K\c-2  V^p-HK^  V^p  =  0. 

The  electromagnetic  term  has  been  dropped,  since  it  is  seen  to 
equal  zero  by  virtue  of  (3),   Next  we  observe  that 
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K^  ^p^a  -   ^a  ^p^""    • 


Hence   the   preceding   expression  may   be   written 


or 

^■^  o  p  ci  pa 

If,  in  addition  to  having  assumed  (3),  we  now  also  use  (l),  it 
follows  from  this  last  expression  that  the  four  equations  repre- 
sented by  (7)  are  linearly  dependent.   Hence  (7B)  is  a  consequence 
of  (7A),  (7C),  (1),  and  (3). 

The  term  in  ( 7A )  which  contains  the  Lagrangian  time  deriva- 
tive is  identified  as  the  time  rate  of  change  of  momentiim  per  unit 
volume  of  the  fluid  element  (i.e.  moraentian  and  volume  seen  by  an 
observer  moving  at  velocity  -  u  with  respect  to  the  fluid).   It 
is  interesting  to  examine  (7A)  for  the  case  u«  c.   The  momentum 
per  unit  voltune  is  seen  to  become  (remembering  the  definition 

i   =  e   +  £-  ): 
o  -   o   p^ 

Po(l^%^^  )  - 
c    "^o 

e 

Since  the  total  mass  density  is  p  ( 1+  — S) ,  the  inclusion  of  the 

c 

p/p   term  in  the  above  expression  might  appear  strange.   The 

physical  reason  that  p/p   actually  does  contribute  to  the  inertia 

arises  from  the  fact  that  there  is  a  supplementary  flow  of  energy 
across  a  fluid  element  due  to  the  work  done  by  the  pressure  of 
the  element  on  its  surroiindings  as  the  element  moves,   (That  is 
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the  neighboring  fluid  does  work  on  the  receding  side  of  the 
element  while  it  receives  work  from  the  advancing  side,  so  that 
energy  must  be  flowing  across  the  eleraento  )   This  flowing  energy- 
has  inertia  which  resists  a  change  of  flow  rateo   An  interesting 
consequence  of  the  fact  that  the  enthalpy  -  and  not  merely  the 
internal  energy  -  contributes  to  the  inertia  is  that  for  suffi- 
ciently high  pressures  a  relativistic  increase  of  inertia  occurs 
even  when  u/c  ->  Oo 
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6,   Determlnacy  of  System 

We  wish  to  compare  the  number  of  independent  eauations  and 
variables  to  indicate  the  determinacy  of  the  system.   (B  and  H 
will  not  be  coixnted  as  different  variables,  as  we  will  implicitly 
assume  the  constitutive  relation  B  =  ij,H.   Similarly  we  implicitly 
ass-urae  D  =  XE.  ) : 

A,  Four  dimensional  formulation 

(1)  Consists  of  1  eauation  and  l\.   unknowns:  K,  ,  Kp,  K^,  K.  , 

(2)  Consists  of  1  eauation  and  2  unknowns:  p,  p  . 

o 

(3)  Consists  of  3  independent  equations  and  adds  6  new 
unknowns  to  those  already  appearing  in  (1)  and  (2): 

^1  *  ^2  *  ^3  *  ^1 '  ^2  *  ^3  ° 
(ij.)  Consists  of  3  independent  equations  and  adds  no  new 

unknowns . 

(6)  Consists   of  1   equation   and   adds   no   new  imknowns. 

(7)  Consists   of   l\.  independent   equations    and   adds   1  new 

■unknown:    i    . 
o 

Thus  we  have  13  equations  in  13  unknowns,  so  that  the 
system  is  determined  (when  boxindary  and  initial  conditions  are 
prescribed) . 

B.  Three  dimensional  formulation 

(2)  Consists  of  1  equation  and  2  unknowns:  p,  p  . 

( 3A )  Consists  of  3  independent  eauations  and  9  \xnknowns : 

E-^ ,  E2  ,  E -J ,  B-|^ ,  B2  ,  B ^ ,  u^ ,  U2 ,  Uo  . 
([(A)  Consists  of  3  independent  equations  and  adds  no 

new  unknowns o 
(6a)  Consists  of  1  equation  and  adds  no  new  unknowns, 
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( 7A )  Consists  of  3  eouations  and  adds  1  new  unknown:  1  , 

o 

(7B)  Consists  of  1  equation  and  adds  no  new  unknowns. 
Again  this  system  is  determined  since  we  have  12  equations 
in  12  -unknowns. 
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7,   The  Stress  Energy  Tensor  (Alternate  Formulation) 

The  following  definition  constitutes  a  combination  of  the 
mechanical  and  electromagnetic  stress  energy  tensors: 

(8)   tP  .  p(l.o-2i^)Ky.p6P-,-\^Br^(^^)-V=B^,5P  . 

It  is  obvious  that  T  „  =  T„ 

ap    pa 

A  generalization  of  (7)  in  that  it  includes  also  electromag- 
netic mass  and  moment"um,  as  well  as  mechanical,  can  be  written  in 
terms  of  this  tensor  in  the  form: 

(7')     VpT^  =  0. 

Equation  (7'),  while  in  itself  not  equivalent  to  (7),  could  never- 
theless be  used  in  its  place  without  altering  the  content  of  the 
original  system  of  equations.   (7')  is  conservation  form,  (7)  is 
not.   The  reason  the  formulation  would  be  equivalent  lies  in  that 
(7)  and  (7')  differ  essentially  only  in  that  the  latter  contains 
additional  information  derived  from  Maxwell's  equations  appearing 
separately  in  the  formulations  anyway.   Eouation  (7')  is  equiva- 
lent to  the  following  two  eauations: 

P^(c  +e  )+p     p  „)       fpo(°  "^^n^"*"?     RR    vv 
1     c  -u         i^c  J        I       "'^  [^^    _, 


2 


2iic^ 


_  (   Po(c2+e^)+p  4-    r2    p2 
(7'B)  V.{  ^-^o s£_  +  B  ^   E 


t       2   2         2ll    ^   2 

'       C  -U  ^     2^LC 
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.2 


p(c  +e  )+p   ^  p. 


Equation  (7 'A)  describes  the  time  rate  of  change  of  the  combined 
mechanical  and  electromagnetic  momentum  at  a  point o   Equation  (7'B) 
is  the  relativlstic  form  of  the  law  of  conservation  of  energyo 
Clearly  the  quantity 

~^ ■*■   ?fr  "■  7-^ 

T    u  ^        2llc 

1  -  -^ 

c 
is  to  be  identified  with  the  total  electro-mechanical  mass  energy 

density  as  seen  by  an  observer  moving  at  velocity  -u  with  respect 

to  the  proper  frame »   Also  the  net  energy  flow  density  seen  by 

this  observer  is  obviously: 

„  Po(^^+e^)-^P  ^  ExB 

1   ^         ^ 
1   -^ 

c 

The  electromagnetic  parts  of  energy  and  power  density  are 

seen  to  be  identical  with  the  classical  form,  while  the  mechanical 
parts  differ  raarkedlyo   In  particular  it  may  be  asked  why  the 
mechanical  energy  includes  the  p  u  c""   term  in  the  numeratoro   The 
reason  is  essentially  that  an  observer  in  going  from  rest  to  some 
finite  velocity  with  respect  to  a  fluid  element  sees  it  undergo 
a  Lorentz  contraction  while  linder  the  pressure  of  the  surrounding 
medium,  which,  to  him  constitutes  work  being  done  on  the  element! 
thus  the  internal  energy  density  seen  by  an  observer  moving  with 

respect  to  the  fluid  is  increased  more  than  merely  by  the  factor 

2  -2  =1 
(1-u  c   )    when  pressxire  is  involved,,   If  u  «  c  and  (  6a  )  is  sub- 
tracted from  (7'B)  the  classical  form  of  energy  conservation  is 
obtainedo 
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It  is  worth  remarking  that  the  above  identification  of 
total  energy  and  power  densities  as  seen  by  an  observer  in 
relative  motion,  while  supplementary  to  the  already  determinate 
formulation,  is  in  an  important  sense,  not  trivial.   For  the 
quantities  p  ,  e  ,  and  1   (=  e   +  pp~  )  entering  the  formulation 
in  no  way  serve  to  disclose  the  corresponding  densities  of  matter. 
Internal  energy,  or  enthalpy  as  seen  by  an  observer  in  relative 
motion.   Thus  if  the  formulation  is  to  be  related  to  quantities 
measurable  by  an  observer  in  an  arbitrary  velocity  frame  these 
supplementary  identifications  must  be  ass\imed. 
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Part  II  -  Small  Amplitude  Hydromagnetlc  Disturbance  Velocities 
1,   Introduction 

The  system  of  equations  formulated  in  Part  I  will  be  shown 
to  be  hyperbolico   The  characteristic  velocities  of  small  ampli- 
tude disturbances  will  be  obtained.   It  will  be  seen  that  unless 
restrictions  are  placed  on  the  mechanical  parameters  describing 
the  medi\mi,  disturbances  propagating  faster  than  the  speed  of 
light  would  be  possibles 

Computations  of  relativistic  hydromagnetlc  shock  (i.e.  finite 
amplitude)  propagation  velocities  have  been  made  for  certain 
particular  cases  of  relative  orientation  between  the  vectors 
involved  (e,go  for  the  magnetic  field  parallel  to  the  shock 
front;  or  for  the  material  velocity  jump  vector  parallel  to  the 
shock  front  normal,  etc„)  by  de  Hoffmann  and  Teller*-   .   In  the 
present  case,  while  treating  only  small  amplitudes,  we  will  con- 
sider characteristic  velocities  of  disturbances  of  all  permissible 
orientations  between  the  field  vectors. 
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2o   Equations  and  Assumptions 

The  non-relativistic  coiinterpart  of  the  results  to  be 
obtained  here  has  appeared  In  Reference  [?]  (loeo  for  the  stan- 
dard Lundqulst  equations).   Since  It  will  be  interesting  to  com- 
pare the  present  results  with  those  of  Reference  [?],  it  is  con- 
venient to  employ  the  three-dimensional  relativistlc  formulation 

of  Part  lo 

The  three  dimensional  relativistlc  formulation  consists  of 
Equations  (2),  ( 3A ) ,  (I4A),  (  6A  ) ,  (7A),  (7B)o   In  addition  to 
these  equations  we  assume  the  followingt   only  small  amplitude 
waves  are  to  be  considered.   We  consider  only  the  proper  frame 
of  the  fluldo   (If  a  distxirbance  advances  more  slowly  than  light 
in  the  proper  frame,  it  will  also  be  slower  in  all  other  Lorentz 
frames. )   The  charge  density  in  the  proper  frame  is  assumed  to 
be  zero. 
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3.   Equation  for  Characteristic  Velocities 

Using  (3A)  we  eliminate  the  variable  E  from  the  other 
equations  of  the  formulation.   Within  the  context  of  obtaining 
characteristic  velocities  ,  we  may  conclude  that  Equation  (2) 
implies  the  relationship 

dp       _  dp    o 
3x, 


"■a 


%  9xa 


Using  this  Vp  is  replaced  by-r^Vp  in  (7A).   Then  the 

Po    ° 
remaining  equations  of  our  formulation  become: 


(I^^)  V^H  +  Vx(H  X  u)=  0 


(6A^) 


V. 


'1  - 


u 


+  V 


=     0 


(7A^) 


f^    V  Po+  liHxfVxH-  Vt  X,tiHxuj  - 


-    (tJ,Hxu)(  V  •  Xt^Hxu)   =  0 


(763^) 


{1+-|) 


1^  +    (ixHxu).    IV^-  Vt^t^H^4=   0- 


""'The  characteristic  velocity  formalism  to  be  employed  is  in 
effect  as  though  the  discontinuities  correspond  simply  to  the 
limiting  case  of  very  large  space  and  time  derivatives.   Thus, 
within  such  a  local  x,y,z,t  region  of  sharp  (but  slight)  change 
of  the  state  variables,  the  "operator  of  the  equation  of  state" 
(equation  2)  is  comparatively  space  and  time  independent. 


-  28  - 


In  the  proper  frame  (u  =  0)  (7B-,  )  becomes,  assuming  small  ampli- 
tude disturbances: 

(7B  )      !io   ^  1_3£ 

This  relation  does  not  contribute  to  the  set  of  simultaneous 
equations  needed  to  determine  the  characteristic  velocities 
(since  a  linear  combination  of  both  time  and  space  derivatives 
is  reauired).   Further  since  we  wish  to  carry  the  1   variable 
along  in  the  formulation,  ( 7Bp )  (which,  in  conjunction  with  (2), 
states  that  1   is  a  function  of  p   only)  will  not  be  used  to 
eliminate  1  .   Thus  we  do  not  use  (78^).   Assuming  small  ampli- 

o  i- 

tudes  and  the  proper  frame  for  Eauatlons  ( IjA,  ) ,  (6a,  ),  ( 7A,  ) , 

the  formalism  of  the  theory  of  characteristics"  leads  to  the 

relations 

(IjAp)      -v5H  +  H5u^  -  H^5u  =  0 

(6A2)      -v5p^  +  p^5u^  =  0 

(7A2)   -p  (1+  -|)v5u+a^n5p^^-Hn(H•5H)-^iH^5H+>C^L^vH(H•5u)-;i;^l^vH^5u=0 

0 

where  v  =  characteristic  velocity  (a  number) 

n  =  xinlt  normal  vector  to  wave  front 

H  =  n.H 
n 

u^=  n.u 

^Po 
5(  )  indicates  jump  (small  In  this  case)  in  a  variable. 

Equations  ( IjAp  ) ,  (  6Ap  ) ,  (7Ap)  parallel  closely  their  non- 
relatlvistlc  counterpart  presented  on  page  10  of  Reference  [7]» 

"See  for  example  Reference  [7]. 
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In  fact  (l+Ap )  and  ( 6Ap )  are  seen  to  be  identical  respectively 
to  (B^)  and  (B^)  of  Reference  [7];  while  (TAg)  differs  from  (B2) 
by  a  factor  in  the  first  term,  and  also  in  that  the  last  two 
terms  in  ( 7Ap )  do  not  appear  in  (Bp)o   The  derivation  of  these 
two  additional  terms  is  given  in  Appendix  1, 

Substituting  into  (7A2)  for  5H  (via  (I+A2))  and  for  5p^ 
(via  (6A2) )  gives 

(7A3)    v^  [[p^(l+  ^)+}^ti.V]5u-H[;^ti2(H.5u)]j 

c 

=  p^a^n5u^+ti,n(H.[H5u^-H^6uj)-^LH^(H5u^-H^5u)  . 

Equation  (7Ao)  constitutes  3  homogeneous  equations  in  the 

variables  5u  ,  5u  ,  5u  .   The  determinant  of  this  set  of  equa- 
x'   y*    z  ^ 

tions  is"  _ 

i       n   '^  o  H, 


-p^v  (1+  -^)+PQa  -^       p^a  -j/H  -H^ 
C         n  Jl 

2   2 


0 


0 


0  '   0        -p  v^(l+  -1+  ^t"  ^  )+^iH^ 

o      ^2    p     ^  n 
c    "^o 

This  is  a  polynomial  cubic  in  (v  ).   (The  determinant  of  equations 
(7Ao)  should  be  essentially  the  same  as  the  determinant  for  the  7 
equations  constituting  (1^2),  {6a^),    ilk^)    -   except  for  the 
possibility  of  zero  roots  (contact  discontinuities).   This  possi- 
bility is  readily  examined  by  putting  v  =  0  into  ( I4A2  ) ,  (  6A2 ) , 
(7Ap),  and  observing  that  the  resulting  relations  are  in  fact 
permissible  for  appropriate  parameter  values.)   Writing  it  out 
"^"Por  a  derivation  of  this  determinant,  see  Appendix  2, 
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and  equating  to  zero  gives  the  equation  for  the  3  values  of  v  , 


-  V 


[(p^a2+tiH2)[l+  -°]+   a.^El)C^^y   ^^E^) 


If  we    set   equal   to    zero   the   relativistically  Introduced   quantities 

— §•  and  X[i  H      (due    to   displacement   current),    ( 7Ai  )    reduces    to    (2,7) 

c  ^ 

on  page  11  of  Reference  [71 •   (Except  that  the  latter  has  an  addl- 

2 
tional  root  of  v  =0.) 
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l+o   Examination  of  Characteristic  Velocities 
Ao   The  system  is  hyperbolic 
The  last  factor  of  the  determinant  (7Ai  )  Is  quadratic  in 
(v  )  and  has  the  discriminant 


•l|ac  E 


a 


+2 


''a2(p^[l+  ^]+eI%/) 


m.h2[i+  -%] 


2  „Z 


-l;(p^[l+  -|]  +  H^)Cti,^)[l+  -iKa-^^H;). 
c  c 

It  is  readily  seen  that  the  discriminant  cannot  be  negative, 

2 

It  is  thus  fiorther  seen  that  both  solutions  for  v   obtained  from 

the  last  factor  of  the  determinant  are  positive  real„   Clearly 
all  6  characteristic  velocities  obtained  from  {7Ai  )  are  realo 

Bo   Intermediate  wave  cannot  exceed  speed  of  light 

2 
The  factor  linear  in  (v  )  of  (2D)  gives  what  corresponds  to 

the  "intermediate  wave"  of  Reference  [']]'. 

r2 


(7A^) 


liH, 


v 


n 


INT. 


2..2 


p^[l+  -4]+>tti.^H 


Because  of  the  additional  terms  in  the  denominator  in  the  pre- 
sent case  (as  compared  to  Reference  [?],  we  see  that  the  inter- 
mediate wave  cannot  exceed  the  speed  of  light.   (Speed  of 

-1/2 
light  =  (  X\i)      '       =  c,  in  the  present  formulationo  ) 

Co   Fast  wave  not  identically  slower  than  light,  but  slow 
wave  is 

2 

To  examine  more  concisely  the  two  solutions  for  (v  )  obtained 

from  the  last  factor  of  the  determinant  we  Introduce  the  notation 
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R  =  p  +  dii_4_  ,    r  =  p  +  J- ,  •  p^  <  r  <  R, 

c  c 

2         2 
Then  the  solutions  for  Vp.gm  and  Vg^^^  are 


(7A,)     v^ 


2 
a 


r+^iH^±  /(a^r+tiH^)^-l4Ra^M-H^ 


PAST  1 

SLOW  2R(1+  -|) 

c 

Prom  (7A/)  It  Is  seen  that  the  slow  wave  velocity  cannot  exceed 

the  speed  of  light o   (To  see  this  compare  Vg^„^  with  the  corres- 

2   2 

ponding  expression  obtained  if  the  term  -l\Ra   [iH   in  the  radicand 

2   2 

is  replaced  by  the  term  of  larger  magnitude  -Ijra  ^.H  .) 

The  fast  wave  velocity  is  not  greater  than  the  speed  of 

-  2  - 1  /2 
light  unless  a(l+i  c   )   '   is  greater  than  the  speed  of  light, 

2 
To  show  this  we  solve  (7A/-)  (upper  sign  option)  for  a  ,   The 

result  is: 

2    2   ,,,  ^o,  ^pastPo^I''  ^^  -^  H(v2^S^/C^2-^,) 


PAST'     2'   p         i     P   P       „ 

^PASTPo(l"-f^<(4sT^^  -^^ 
c 

This  expression  immediately  indicates  that  if  v„.„_  >  

r  Ao 1     ^ • 

a         1  ^^^ 

then    ^     >  -^  , 

c 

Special  Cases  of  (7A/) 

Some  of  these  will  illustrate  that  the  speed  of  light  is 

2 
exceeded  if  a   is  sufficiently  large: 

If  H      =  0  o 

n  ^        tt2 

2  ^  Po  2  ^   ^ 

^PAST  i  ..    2^2  ^SLOW 

-,      _o        /i[i  H 

2  p 

c  ^o 
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If  H   =  IhI 

S -1  2 

uH 


2„„     =^     v^        -      '° 


SLOW  i       "PAST  1     2^2 

(OR  PAST)    1+  -|     (OR  SLOW)    1+  _|+ MJi_ 
c  c     *^o 

The  choice  of  upper  or  lower  Index  in  the  last  two  expressions 
(indices  coupled)  depends  upon  the  comparative  numerical  values 
of  the  respective  expressionSo   The  equation  of  state  of  the  gas 


must  satisfy^  do  ~  ^  "^  ^^^   "^  ^o^^   "   "'"^  ^'^®  last  expression 


(Alfv^n  wave  case)  H   =  |HJ  to  first  order  only^ 


3k 


APPENDIX  I 

Evaluation  of  Expression  Corresponding  to 

Bx— (Bxu)  for  Characteristic  Velocity  Formalism 
dt 

Bx^(Bxu)  =  i(By  ^(Bxu)^-  B^  ^(Bxu)y) 

-    J(  ) 

+  k(  ) 

+  j(  ) 

+  k(  ) 

8u 


Assioming  small  amplitudes  we  have  ^Bj^'^^  =  B^  aT^*   That  is 
the  B  components  here  act  as  constants  with  respect  to  the  time 
derivative  operator.   Using  this  fact,  we  get  after  some  manipulation; 

BX  I^(BXU)  =   KByB^  ^   +  B^B^  _^  .  b2  5^  +  32  ^) 


A 


+  k(  ) 

Bx^(Bxu)=   5(B^(B.  fH)  -g^B^) 

^  3'«y'^-  If'  -  ^  ^'' 

^  i^'«z(B-  I?)  -  ^  3') 

Bx^(Bxu)  =B(B.  if)  -^''i- 

Under  the  formalism  for  characteristics  vje  replace  qt-  by  -v5u. 
Hence  Bx-^(Bxu)  corresponds  to:  -vB(B»5u)  +  vB  5u, 
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APPENDIX  II 
Obtaining  the  Determinant  of  the  Three  Equation 
System  Constituting  the  Relation  (7Ao) 

Choosing  the  field  H  to  be  in  the  x-direction  (to  first 
order)  the  three  equations  comprising  the  relation  iJA-.)   may  be 
written  (after  dropping  higher  order  terms); 

X  component   p^v^(l+  -|  )5u^  =  p^Su^la  nl+p^6v.^j&   n^nyj"^ 

c 

y  component   p^v  ( 1+  -|  +  -^ — )5Uy  =  Po5u^|a  nyn^j+p^5Uy|a  n^ 

2  2 

+   — -in'^+n    )  r  +p  5u  .^a  n  n  +  — -   n  n  I 

1      2  2 
z  component   P^v^d^  -°  +  >^ii-H_)g^^  ^  Po^^xP^z^xj-^-Po^^yK^z^; 


2         ^  tt2 

+ 


^}^Po-.{^^"^%<"^"^} 


We  chose  H  to  be  in  the  x-directlon.   We  lose  no  essential  generality 

if  we  further  choose  n  =  Go   (Since  the  wave  front  normal  can  still 

z 

be  at  an  arbitrary  angle  to  Ho)   Putting  n  =  0  in  the  above  three 
equations,  the  determinant  becomes 

-PoV^(l+-f)-^  Po^^^x         Po^^^x^y  ° 

c 

Po-Vy  -Po-'(l^  7  -^  ^Hp^.^W  0 

,0  0        -p^v2(l+  ^  +  ZJ^)+^^nl 

c      o 

H 
Since   H    is    in   the   x-direction  we    have   n     =  tttj  <>     Also,    since   n     =   0, 
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h2 

we  have  n^  =  1-  n^  =  1-  "4.  Expressing  the  n^  and  n  in  the 

7       -^      H 
above  determinant  in  terms  of  H  and  H^,  it  becomes  the  determinant 

given  in  (7A.  ) . 
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